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A new way to write the massive scalar and fermion propagators on a background of a 
weak gauge field is presented. They are written in a form that is manifestly gauge-covariant 
up to several additional terms that can be written as boundary terms in momentum space. 
These additional terms violate Ward-Takahashi identities and need to be renormalized by 
appropriate counterterms if the complete theory is to be gauge-covariant. This form 
makes it possible to calculate many amplitudes in a manifestly gauge-covariant way (at 
the same time reducing the number of Feynman diagrams). It also allows to express some 
counterterms in a way independent of the regularization scheme and provides an easy way 
to derive the anomalous term affecting the chiral current conservation. 
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1 Introduction 


In a quantum field theory, when one wants to calculate a specihc amplitude 
in the form of a formal series, one often uses Feynman diagrams [1]. This 
is doubtlessly an amazing method which allows to swiftly write a number of 
expressions contributing to the desired amplitude. However it is not with¬ 
out weaknesses. Because it is based on a division of the lagrangian into the 
free-held lagrangian, the interaction lagrangian and the counterterms, the 
individual expressions often do not have the symmetries of the theory if the 
parts of the lagrangian do not posses them on their own. One such symmetry 
that is never manifestly preserved is the gauge symmetry, since the initially 
gauge-invariant terms in the lagrangian are broken up into pieces, and differ¬ 
ent pieces are treated differently. Moreover, if we consider the gauge held to 
be not only the classic, background held, but a dynamic quantum held, we 
need to add the gauge hxing term to the lagrangian so that we could con¬ 
struct a propagator for the gauge held. Eventually the individual expressions 
are not gauge-covariant, and even when summed up the gauge-covariance of 
the hnal result is often not visible at once, and is only recovered because 
of cancellations between many terms. Sometimes even the Green functions 
themselves aren’t gauge invariant, and only the physical observables, like 
cross-sections and scatering amplitudes, are. 

The alternative approach by Schwinger [2] and Tomonaga [3] is gauge 
covariant, but at the same time it is difficult to calculate with and opaque. 

While the development of the numerical method in last decades and rising 
power of the computers made it possible to calculate many diagrams with 
relative ease, at the same time, by receiving only the hnal result we lose the 
insight on how it came to be. Algorithms that perform the calculations rarely 
even give the result in the form of a gauge-covariant formula, producing in¬ 
stead only numerical values that one needs a lot of experience to comprehend 
easily. I believe that to preserve the understanding of the gauge theory one 
should look for analytical expression if possible, and be more aware of how 
certain results are produced. However, the method I’m going to present is 
also algorithmic and can be put on a computer if you wish so. 

I am going to present a way of obtaining the formulae for scalar and 
fermion propagators that are gauge invariant and can be used as building 
blocks in Feynman diagrams. It bases on the background held method HEIE] 
used to calculate the ehective action. That method already guarantees that 
the hnal result is gauge-covariant, but in this work I want to formulate such 
Feynman rules that not only the hnal result, but also intermediate steps are 
given by gauge-covariant expressions. It is done by expressing the propaga¬ 
tors of the particles on the background of a gauge held using path-ordered ex- 
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ponentials and a quasi-local formal series made of exclusively gauge-covariant 
quantities. One may be familiar with the following approximation for the 
quark propagator in a background of gluon fields: 


{q{xi)q{x 2 )) ^ Pexp ( - 0 / G“(x)T“ck^) {q{xi)q{x 2 ))o ( 1 ) 

J X2 

where Pexp denotes a path-ordered exponential (along some path connecting 
points Xi and X 2 ) and subscript 0 refers to free quark helds. Typically, one 
needs to average over all possible paths, the concept that is used in the lattice 
gauge theory introduced by Wilson [7]. The idea I follow in this work is to 
focus one path - the straight line - and hnd the corrections to the formula 
above in the form of a series constructed with gauge-covariant operators. 

First of all, the corrections I hnd exhibit nontrivial spinor structure, the 
true propagator cannot be expressed simply by multiplying the free propaga¬ 
tor with some function of the gauge held. Notably, amongst other terms in 
the propagator, I hnd one proportional to 75 , which can be used to obtain the 
chiral anomaly (ABJ-anomaly [HIIS]) in a new, transparent way, without the 
need of any regularization. The full formula for the propagator, in the form 
of a series, can also be used in Feynman diagrams, making them to produce 
only gauge-invariant expressions and sparing us the need of calculating terms 
that would cancel each other in the hnal result anyway. 

In this work, I construct the formulae for the scalar and spinor propaga¬ 
tors. For the full set of rules, we would also need the rules for the lines of 
gauge bosons (which in the backgound held method never appear as external 
lines, but can appear inside the diagrams). However, this case is much more 
complicated, especially if we want to consider a general gauge (which may 
also include ghosts) or spontanous symmetry breaking like in the SM. Since 
even with only scalars and fermions we can see some interesting features of 
the method I present, and because of the limited space I have, I decided to 
leave the vector boson case out of this paper, and I plan to publish it in a 
following work. 

2 Presentation of the method, scalar case 

Let us consider a theory of a scalar held minimally coupled to an external 
gauge held: 

C = D^^(t)^D^(t) - mVV (2) 

where 

Df^cf) = - A^(j) 
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A„ = -igAir = -A* ( 3 ) 

For now, we are going to consider the field to be a background field, so we 
do not include its kinetic energy in the lagrangian, and it will not appear as 
a propagating field inside the diagrams. This lagrangian is invariant under a 
local redefinition of fields with U{x) = exp(—i5f0“(x)T“): 


(j){x) = U{x) 4 >'{x), Afj,{x) = U{x)A'^{x)U{x) ^ + d^U{x)U{x) ^ ( 4 ) 

In the path integral formulation of the quantum theory, we have 


Let us consider the 2 -point Green function 


J=0 


This function should follow the appropriate gauge covariance rule 
G{xi,X2)[A] = U{xi)G{xi, X2 )[A']U{x2)~^ 


( 6 ) 

( 7 ) 


which gives us the following Ward-Takahashi identity: 

d SG^{Xi,X2)[A] rabc Ab( JG^{XuX2)[A] ^ 

dxf^ 6 A-{x) ’ dA-^ix) 

= -ig (G^{xi, X2)[A]T''6 {x2 - x) - 6{xi - x)T‘^G^{xi, X2)[A]^ ( 8 ) 


The issue we are going to address is the fact that when we calculate Green 
function G2{xi,X2)[-A\ using perturbative expansion, the satisfaction of eq. 
(|^, which comes from the gauge-covariance, is completely invisible. The 
formulae we get are not-gauge-covariant, only when summed up the gauge- 
covariance is restored. As a consequence, the Feynman diagrams that use 
free propagators of 0 and add interaction with gauge field as a perturbation, 
don’t give gauge-covariant results as well, unless summed up. This means 
that during the calculation of every single Feynman diagram, we needlessly 
calculate also some not-gauge-covariant expression that is bound to cancel 
out with other such expressions from other diagrams. In this work we are 
going to show that it is possible to present function G2{xi,X2)[A\ in gauge- 
covariant form, and using that form to simplify some Feynman diagrams. 

The method focuses on using the function G{xi,X2)[A\ in the form 

G^{xi,X 2)[A] =Vexp (^J A^G^{z-,xi,X2)Vexp (^J A^ ( 9 ) 
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where Pexp(f A) denotes a path-ordered exponential of the gauge held along 
a straight line, which can be dehned by the following differential equation; 


^"Pexp = a>^A^{x Xa)Vexp 


Pexp /;A =1 


( 10 ) 


It can be shown that it satisfies 

f*Xi 


Vexp ( / A^ = U{xi)Vexp ( / A'']u{x2)~^ ( 11 ) 

^ J X2 ' ^ J X2 ' 


or 


d '^^exp X 


r'Xl 

fX 2 


gf^Aiix)- 


iPexp f 


pXl 

iX 2 


5 A^^{x) 


dxf^ 6A^{x) 

= -.g(v^,{rA)T‘5(x.-x)-5ix,-x)T‘Ve.,{ 


f*Xi 


X 2 


( 12 ) 


Combining equations (^ 8 |-( 11 ) we get 


G{z-,Xi,X2)[A] = U{z)G{z-,Xi,X2)[A']U{z) ^ 


( 13 ) 


d 6 G^{z-Xi,X 2)[A] _ rabc Ab( JG 4 ,{z]Xi,X 2 )[A] ^ 
dxf^ 6A-{x) ’ M^(x) (14) 

- i^f [G^{z; xi, X2)[A],T'^]6 {z - x) 

The simplihcation we have obtained with respect to eq. ([^ gives us hope 
that G(j){z; xi,X2) can be written in simpler form than G{xi,X2). Especially 
if we focus on weak external helds A, and calculate G^{z;xi,X2) as a series 
in powers of A{z) and its derivatives, we expect to see only gauge covariant 
structures, like A^^,, D^A^p etc. 

The point can be in principle arbitrary, but the easiest choice seems to 
he z = Xi, z = X2 oi z = though if we know beforehand in what 

Feynman diagram they will be needed, a better choice may be available. 
We will present the results for z = as in general, it is the one which 

simplihes many Feynman diagrams the most. In this case I’m going to write 
the function Xi,X2) as a function of two arguments Xi — 

X2), from now on. 
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To calculate G^{x, a), we move the path-ordered exponentials from ([9j) to 
the other side: 


G^{x, a) [A] = Vexp 


X + -a, X — -a)Vexp 


A) ( 15 ) 


We will use G^{x + ^a,x — |a) in the form of series in A, which comes 
naturally if we use Feynman diagrams. At the leading order we calculate 
only tree-level diagrams, with one scalar line and an arbitrary number of 
gauge field lines attached to it. 



Fig. 1. The tree-level diagrams that contribute to the scalar propagator in an external gauge field. 


We also need to expand the formula in powers of moment a/derivatives of field 
A. Let us remember that this is just G(a;i,X2), which doesn’t yet transform 
simply under gauge transformations, and without path-ordered exponentials 
we don’t expect to be able to gather different terms to create only gauge- 
covariant tensors like A^^, = S^Aj, — — [A^,Ai,]. For now, let us just 

list all the terms that appear up to the order of ( 9 (A^), (counting derivatives 
and A to be of the same order, as suggested by the form of the covariant 
derivative 


G^{x + a/ 2 , a; - a/ 2 ) [A] = 

^ ^ (p2_^2)2^''(^) + 


9 piv 


(p2 


m 


2\2 


+ 


-^P^lPv 


(p2 


m 


2'|3 


''^9apP^l , -2ipaPf)Pi, 


+ 


(p2 _ 77^2^3 (^p2 _ 


M{x)K''{x)+ 

d^d^M{x)+ 


V 

Jo 


r 

J {2tt)0 
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1 [yA>-(x),A-(x)] + 


(p2 _ Jji2'j3 (jp. _ jyPy 

[+ Ppg^,v) ^ -^iPpPuPp 


\ {p^ — m?Y 


(jP _ jypy 


M{x)hy{x)K^{x) + 


29a09pi/ PaPpgpv PpPa9ul3 PpPl39va ^PpPu9o!fi 

I / O O \ /I I 


— m2)3 (p 2 _ 77^2^4 

+ P^){3“9^^‘‘(x),A‘'(x)}+ 

9af39pi' Qi'agp^ ^ ‘^PaPp9pi' ‘^PpPjSPi'ct ‘^PuPaQpp ‘^PpPu9afi 


(jP _ jpPY 
X (a“A^(a;))(a^A"(a;)) + 


(jP _ jypy 


X 


9pa9pv -‘^PpPa9pu - ‘2ppPa9up “ ^PpPu9pa “ ^PpPp9v 


— 77^2^ 


— mP'Y 

+ (9"a'‘(x))a''(x)a<'(x)+ 

9pa9pu 9pa9up ^ ‘^PpPa9up ‘^PpPa9pu ^PpPu9pa A ^Pi/Pp9po 


X A^(a;)(a"A'^(a;))A^(a;) + 




X 


9pa9yp ‘^PpPaQup + ‘2.PpPa9pv + ^PpPp9uci + ^PuPp9pa 

I / O r» \ /I I 


— m 2 p 

'^GPpPuPuP, 


+ 




9py9po 


+ 


(p 2 _ 7772^3 

(p 2 _ 7772)5 

(9(((9, A)^) + (boundary terms) 


— w?Y 

K^^{x)^y{x)^'^M{x) + 

-^{PpPaPpy + PpPaPyp + 

(p 2 _ 77^2)4 

A^(a;)A'^(a;)A^(a;)A'^(a;) 


(16) 


The boundary terms in this expression come from the ambiguity of redehning 
momenta in the derivation of this result; were all the integrals convergent, 
one could make the redehnition of integration variables like 


Pi^P - 
P2^ P+^q + Xq 


(17) 


for any value of A, and obtain the same result. However, since some of the 
integrals are not convergent, we receive several (hnite number) terms of the 
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form 


(18) 


A,; 


A 

dpu 




where A* are, at this point, arbitrary constants. We will get even more bound¬ 
ary terms in the formula for G^{x, a) because the path-ordered exponentials 
from (15) produce many terms that will need to be integrated by parts. For 
any operator O we have 


Vexp ^ a) OPexp J " A j = 

= O - ^a/,{A^(a;),0} + - [(?^A''(a;), O] + {A^(a;), {A‘'(a;), O}} j + 

+ - {a^a"A^(a;), O} + [a'^A"(a;), {A^(a;),0}] + 

+ 2{A'^(a:), [a'^A^(a:), O] } - {A'^(x), {A‘^(a:), {A^ O} } }) + 

+ ^a^a,apa^(^ - O] + {d^d^AGx), {A'^(a:), O}} + 

+ 3 { A'^ (x), { (x), O } } + 3 A^ (x), [a^A" (x), O] ] + 

- [a^A^(a:),{A^(a:),{A'^,0}}] - 2{A'^(a;), [a^A'’(a;), (A*", O}] } + 
-3{A^(a;),{A'^(a;), [a^A‘^,0]}} + 

+ {A^(a;), {A"(a;), {A^ {A^O}}}}) + ... 

(19) 


Then we use integration by parts to get rid of factors; 

ttp j e“‘^“/(p) = J e“'^“ boundary term) (20) 

As it turns out, with a general regularization|^ these boundary terms cannot 
be neglected and play a crucial role in some diagrams. They are also related 
to certain counterterms that appear in the process of renormalization, an 
example is presented in a section 4. below. For the list of boundary terms, 
see the appendix. 

Eventually, many terms in the formula for Grf){x, a) turn up to be vanish¬ 
ing, and the remaining ones can be grouped together to form gauge-covariant 

^In regularizations that do not violate Ward-Takahashi identities, like the dimensional 
regularization, these terms will vanish and can be skipped early. In general case though 
they need to be remembered. 
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quantities, as expected: 


G^{x,a)[A\=i 


p2 _ ^2 




lipu 


(p 2 _ 17^2^3 


1 / 

-{A^,(a;),A,.(x)} (^(^2_^2)3 “ (p2 _ ^2)4 


+ 


+ (terms of higher order in held A or its derivatives) 
(boundary terms) 


+ 


( 21 ) 


Let us behold how much simpler it has become, compared to (16). The 


reason for this is that there aren’t many gauge-covariant structures that can 
be written at this order, and some of them (like A^^,) cannot couple to a 
scalar because of the lack of proper Lorentz-invariant structures. 


3 Fermionic case 

The same thing can be done for fermions. We start with 

C^[A] = i'lpYD^ip - m'lljil) 

Z^[7],ri,A]=A^AvX,A] = f ^^^ifd^C^[A]+vP+^r,) 


GA^i,X2)[A] = 


. 6'‘W[ri, fj. A] 


i'i](xi)iri(x2) 


( 22 ) 

(23) 

(24) 


i ]=0 


Using the same method (we shall skip the intermediate steps), we can hnd 
that 


G^{x + ^,x-^)[A\ =Vexp 


with 

G,p{x,a) = 


AjG^pix, a)Vexp 


(25) 


= 1 / e 


—ipa 


+ m 


pZ _ 


2 i 

—]D^Ai,p{x) 


c^Apui^x) 


7 ^ Pq, -|- 7 ^ m 


i 2 — m 2 2 (jp, — 


ixp 


ipl-^lpp 


\(jA — m?Y {p^ — m?)^ 


7 


glJ-Ppt' 


(jA — mAY 


-|- m) 


+ 
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1 


9‘ 




Ap^p'-' 


(p 2 _ ^2)4 


+ 7^'"m) + 


^ ^ (7"">a + r^rn) - , / + Y^rn) + 


(p 2 _ 77^2^4 
rpP/y^P^ _j_ p^/yPP^ 


(p 2 _ 7772^4 


(p2 — 777,2)3 

“ 1 “ ■[ (x), Apg-(x) j" 

gt^Pp’^p^ 

- (p 2 _ 77^2)4 

A [A^jj/( x), Ap(j (x)j 


+ 


\ylii/paap^ 7. gPp^pU^a _j_ pCr^!^) 


8 (p 2 — 7772)3 


+ m) 


+ 




+ 


pPpP 


(p2 _ 777 , 2^3 ^p2 _ 7772^4 


2 (p 2 — 7772)3 


(7"">a + 7""^) + 


+ (terms of higher order in field A or its derivatives) j - 
(boundary terms) 


(26) 


where 7 ^^ = 7 ^^^ = etc. This formula is more complicated 


than ( 21 ), because now all possible gauge-covariant structures appear, but 


still much simpler than the fermionic analog of (16) is. 


We would like to stress that this formula is valid in any dimension, and 
it’s independent of the definition of 75 . The only relation between gamma 
matrices that is necessary to derive it is their anticommutation relation, 
{ 7 ^, 7 ^} = 2g^^l. If the dimension is given, it can be simplihed because of 
the fact that sufficiently long antisymmetrized products of gamma matrices 
are 0, for example, in 4 dimensions we hav^y^^^”"" = 0 and 'yp'^p^ = ie^^^^'ys. 


4 Boundary terms and counterterms 


As mentioned before, the boundary terms that appear in ( 21 ) and (26) can¬ 


not be neglected in a general case. Most of them are related to the coun¬ 
terterms necessary to make the effective theory of gauge field finite and/or 
gauge invariant. For the example, let us focus on the fermion theory in four 
dimensions. 


From (23), using explicit form of the lagrangian given by the (22), we can 


3 with the convention that 75 = 17^7^7273^ £9123 = —e °323 _ 
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derive the equation 

SW^[ri,fi, A] 


SA^(x) 


Tj=0 


9Tl{7'‘r“(3^(i,o = 0)| 


(27) 


One of the boundary terms that appear in (26) from integration by parts 

d 


is 


Gp(a;,a) D 


dpp 


^-ipa_ 


f + m 

p2 _ ^2 


We can see that this term will give a contribution to (27): 

6W^[r],fi,A] o X f d 


5Al{x) 


D -2ig^Al{x) / — 


)?=0 


ph- 


dpi, _ 77^2 


which means that 


^-SMl / A;(x)A°‘‘(x) 

rj—0 Z J 


with 


6Ml = -i^ 


d 


pH 


dpP _ 777,2 


(28) 


(29) 


(30) 


(31) 


To maintain the gauge invariance of the effective theory of the gauge held, 
we need to add a counterterm to the lagrangian (22): 


^ + c 

ct 

C,t 3 24“71“^ (32) 

It is exactly what we could get calculating the radiative correction to the 
gauge held mass from 1-loop diagram (Fig. 2.). However, the method pre¬ 
sented in this paper shows that the value of the counterterm is already con¬ 
tained within the coefficient to a boundary term in fermion propagator cal¬ 
culated on the tree-level. 



Fig. 2. A 1-loop diagram that contains a contribution to the gauge field mass term. 

The gauge symmetry will be preserved and there will be no need for such 
counterterm if the boundary terms vanish in a given regularization, for ex¬ 
ample in dimensional regularization. However, other regularizations can give 
non-zero results. 
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5 Boundary terms and axial anomaly 

In this section we are going to use the results obtained to show the origin 
of the ABJ-anonialy [HI [9] in 4 dimensions. For simplicity’s sake, we shall 
restrict ourselves to the case of abelian gauge group but the 

calculations can be performed in a general case. With the propagator in the 
form (26) it is easy to calculate 

iM^)) = = 

r^{x,x) \ = -Tr 


= -Tr 


{x, a = 


(33) 


Up to the terms of order 0{{d, A)^), the only term from G^{x,a) that can 
contribute to this trace is 


G^{x,a) D m F^G^)Fp„{x)y 


fii/pa 


^-ipa_ 


— 1 


4(p2 — 


Assuming that 75 is dehnecQ such that 

Tr{7^^^'"75} = 4ie^^^‘" 


(34) 


(36) 


we hnd the result for (33) to be 

iMx)) = 

= -e"^ m FpGGFp^i^) Tr{ 757 '"‘'^'"} ■ 


-1 


1287r2m2 


+ omAf)^ 


(36) 


le 


e>^’'P^FpGx)Fp^{x) + OGd,Af) 


327r2m 

Independently, we can calculate 

= -Tr| 7 '" 75 G',/,(a:,a:)| = -Tr| 7 ^ 75 G^(x, a = 0 )| 

This time the only terms that can contribute are some of the boundary terms 
G*(i, a) D e" A^(x)a,A,(x) + (62 - x 


(37) 


X 


A. 

dpo 


^-ipa_ 


PP 


{p‘^ — rn?Y 


(38) 


Which is how it should be defined in dimensional regularization, according to’t Hooft 
and Veltman HD], where is defined in such a way that it s equal to 0 if any of its 

indices is different than 0, 1, 2 or 3. 
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where A 2 and A 3 are arbitrary constants associated with the fact that shifting 
the integration variable (momentum p in this case) in an integral that is not 
convergent produces a boundary term, as mentioned before (17), see also the 
appendix. We get 


(39) 

However {J^{x)) needs to be gauge invariant, and to ensure that we must 
either put A 2 — A 3 = — 1 or choose a regularization scheme in which this 
boundary term is equal to 0. Either way 


{jax))^o((a,Af) 


(40) 


Let’s now check the deviation from the naive axial current conservation 
equation: 

{d^J^{x)) - 2im{J5{x)) = 0 (41) 

We can see it’s not satished; instead 

{a.J^ix)) - 2im{Mx)) = + 0((S.A)=) (42) 

which is the well-known ABJ anomaly 0 . 0 . The calculation here was 
made only up to the terms of order of 0{{d, A)^), but from the general 
theory m we know there cannot be any terms of higher order in helds. The 
anomaly in this expression could be avoided but at the cost of losing the 
gauge invariance. It is important to understand that in this derivation the 
anomaly does not come from the divergence of the current, but from (J^lx)), 
which is proportional to 1 /m. 


6 Examples of the application of the method in 
Feynman diagrams 

Expressing the propagators in this form simplifies the calculation of certain 
Feynman diagrams. The most important limitation though is that we’re us¬ 
ing an expansion of G{x, a) in the number of external gauge held lines and 
their momenta, so this method can be used only if the momenta of the gauge 
bosons are much smaller than the masses of particles they couple to. 
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A simple example would be the decay of a light neutral scalar coupled to 
heavy charged fermions: 



(43) 


int 


—Xip'ip'ijj for scalar 

foi' pseudoscalar 


(44) 


Normally, the calculation of the amplitude of (p —)• AA decay would require 
two triangle diagrams at 1-loop level (Fig. 3.), and neither of them is gauge- 
covariant (only their sum is), and the part cancelled by the counterterm is 
hidden within the expression. However, with the method we present in this 
work, we only need one tadpole diagram, which is gauge-covariant by itself, 
with the exception of contribution cancelled by a counterterm which is clearly 
visible (Fig. 4.). 



Fig. 3. Standard 1-loop diagrams that describe tsp AA decay. The propagators of i/j are the free 

propagators in vacuum. 



Fig. 4- The only 1-loop diagram that describe gx A A decay with the method presented in this work. 
The propagator of ijx is the propagator in external gauge field. 


At the 1-loop level, in the case of the coupling without 75 , we have 


out (AA|<p)iji — 
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—8iA 


[ (^^|Tr{A^i^(a;)Apa(a:)}|0)(0|v7(x)|v7) + ■■■ + 

{P Jx 

(p^) l{AA\Tr{A,(x)K{^)}\0} ( 0 |<,(.)|<,) = 

^ ~^ITr{(x)(x)}10)(01V?(x)IV?) + •••+ (45) 

+ (a part cancelled by a counterterm) 


where dots denote the terms with higher number of derivatives of A^. 
The case of the coupling with 75 is even simpler; 


out 19^) in 

= -A /{AA|Tr{75G'^(a;,a;)}|0)(0|v3(a;)|(p) = 

J X 

= -A /"{AA|Tr{75G^(a;,0)}|0)(0|(p(x)|(p) = 


= -ATr{757"^^'^} 




■X 


(46) 


A 


IbTT^m^ 


jp - mlf 

X f {AA|Tr{A^j.(a:)Ap^(x)}|0)(0|(p(x)|(p) H-= 

J X 

(AA|Tr{A'"''(a:)A^^(a:)}|0)(0|(p(a:)|(p) + ... 


Another calculation in this model that becomes much simpler if we use 
formulae derived in this work, is the correction to the mass of (p due to 
background gauge held that appears in the effective action after integrating 
out the '0 field. With this method, it is all contained within a single diagram 
(Fig. 5.) 



Fig. 5. A 1-loop diagram that contains the corrections to the mass term of (f. 


The contribution to the effective action given by this diagram is (assuming 
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the coupling without 75 ): 


‘^eff Z) 2 
1 



^{x)^{y)^i{G^{x,y)G^{y,x)] = 


X Jy 



=-^iA^ / j yj{x + a/2)yj{x - a/2)Ti[G^{x,a)G^{x,-a)] = (47) 


= -2^^' 


(p{x) 



X 


,a)G^{x, -a)} + ... 


For a 7 ^ 0 the boundary terms disappear because of the oscillating factor 
e-ipa^ so the only relevant terms that contribute to this integral are covariant 


terms from (26). The first non-zero terms after the held-independent term 


that we get from standard Feynman diagram are proportional to K^y{x)Kpc{x) 



+ ... = 


+ rn^) 

(p2 - m2 )4 

= (a constant that is a subject to renormalization)-|- 

Therefore the renormalized contribution to the effective action is; 


(48) 


5, 


A 


'eff 3 / 43^2^2 Tr{A^^(x)A^^(a:)} + 


(49) 


With the propagator already in the form (26) the computation is much sim¬ 


pler than it would be to do from the scratch. 


7 Complex diagrams 

In the diagrams in the previous section, in which there were no more than 
two vertices (not counting the interaction with the gauge held), and none of 
them contained the gauge group generators, the formulae drastically simpli- 
hed because the path exponents cancel each other and disappear from the 
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expression. In more complicated diagrams they may remain but even then 
there is a way to combine them to quasi-local, gauge covariant expressions if 
necessary. 


I 

I 

I 



Fig. 6. The simplest case in which we need to remember about Wilson lines. 


If we have more than two vertices, like in the Fig. 6 ., the hrst thing that we 
may notice is that their middle points do not coincide, for example one of 
them can be 

^ 1+^2 

^ ^^ ’^1 - ^ 2 )'^exp J ( 50 ) 

and the other 

/•X2 I ^ ^2+^3 

Pexp A)d( ^^ 2 ~ X 3 )Vexp j A^ (51) 

If Xi 7 ^ X 3 then one of them is expressed in terms of the gauge held in point 
and the other in point To solve this issue, we can use the analog 

of Taylor series in the space with a gauge connection; 

0{x) = Vexp ^ y A^ X ^C>(xo) + (x - xo)^D^O{xo) + 

+ x-o)“(x - Xo)‘'r>„r>„0(xo) + ...) X Pexp (^J a) 

(52) 



We can choose the point Xq arbitrarily, usually one of the vertices or the 
diagram’s "mass center" is a choice that leads to simple expressions later. 
We can also use this formula to "shift" any vertex that contain group indices 
to the point xq. 

This way, if the initial formula contains expressions like 


.X 1 +X 2 . 

Oi 2 { ——)F’exp 


. ^1+^2 
2 


'X2 


A)Id2(a:2)^exp 


f*X2 


f X 2 +x^ 
2 


, X 2 + X 3 , 

023(^y^)... 

(53) 
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they can be written in the following form; 


rci +XO 

S ^ f 2 

.. .O'i2{xo) xVeyiY)[ / AlPexp 

*£l+£2 / 

2 


• X2 


f*X2 


X V^^xq) X Pexp 


f3:o 


UPexp 


'■3^2 . 

AjPexp 

^2+^3 / 

2 


"Pexp ( / A I X 

^ JxQ 

A) xO'3(a;o).. 


'Xo 


(54) 


with O' and V derived from O and V as eq. (52) dictates. As we can see, the 


path-ordered exponentials form triangles. These triangles can be expressed 
in the form of a qnasi-local, gange-covariant series at the point Xq\ 


"Pexp 




'Xl 


Pexp 


'X2 


Pexp 


'XQ 


= 1 - - (^(Xi - XqY{x2 - XqY - {X2 - XqYYi - XQY^^^lu{Xo) + 

- ^( 2:1 +X2- 2xoY{{xi - XoY{x2 - XqY “ (^2 “ XqY X 

+ ... 

Another approach, that wonld let ns avoid such triangles, would be to use 
a general formula with point z = Xq being the same for all propagators 


in the diagram, instead of (50) and (51). Then only the vertices need to be 


shifted with (52), and the necessary path-ordered exponentials are already 


there. However, general formula for G(xo; xi, X 2 ) is more complicated than 
in the case of xq = and contains vectors (xq — X\Y and {xq — X 2 Y- 

Either way, we obtain the formula in the form of a series in the powers 
of the held K^^ixo) and its derivatives. The formula however also contains a 
number of vectors [xi — XjY- If we want to get rid of them to perform the 
integrations over Xi and remain with a single spatial integral, we can chose 
Xq to be some linear combination of xi,.. .x„. Then all these vectors can 
hrst be decomposed into the vectors related to particular propagators, and 


then be turned into the derivatives over momenta using eq. (20) 


While this procedure certainly looks complicated, it needs to be said that 
it is necessary only in the case of very complicated diagrams, with at least 
3 vertices other than coupling to the external gauge held. And even then, if 
we are interested only in the leading contribution for weak held A^j^, or we 
are working with the abelian case, there is a good chance we would be able 
to perform additional simplihcations and reduce the number of the triangles 
of path-ordered exponentials before employing the formula (p^. We should 
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also remember that calculating such diagrams in the standard way is usually 
even more complicated, and one needs to calculate more diagrams to get 
a gauge-covariant result. With this method we work with gauge-covariant 
quantities all long, and we avoid computing some irrelevant contribution that 
cancel between standard diagrams at the end. 


8 Summary 


In this work i present an alternative method of calculating some amplitudes in 
QFT. While it does not make it possible to calculate anything that couldn’t 
be calculated before, it can make some calculations faster, and most impor¬ 
tantly, we can have better control on what happens in intermediate stages of 
computations, since the physically relevant, gauge-invariant terms are clearly 
visible and not masked by the multiple other terms cancelling in the hnal re¬ 
sult. I only calculate what is going to remain, and do not need to consider 
irrelevant terms. Without the need to calculate all the diagrams to get a 
gauge-covariant results, it may also be easier to focus only on some subset 
of them, in situations when we are able to argue that others do not produce 
the contributions that are relevant for the case in hand or that they give 
contributions that are negligible. 

The method has its limitations of course. Because I only managed to 
derive the Feynman propagators in the form of a series in gauge helds and 
its deivatives, it’s only apllicable when the external gauge helds is slowly 
changing, or the momenta of the external gauge bosons are much smaller 
than the masses of particles they couple to. For this reason, some processes, 
like the production of gauge bosons in high-energy collisions, can’t be calcu¬ 
lated with this method unless a better formula for the gauge-held-dependent 
propagators were to be derived. Still, it can be used in many other situations 
and I believe the clarity and simplihcations brought by the application of the 
propagators in the form presented in this work makes the effort of deriving 
them worthwhile. 


Acknowledgements: I would like to thank Adrian Lewandowski, Krzysztof 
A. Meissner and Hermann Nicolai for the discussions and the help in prepar¬ 
ing this paper. 
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A Boundary terms 

For the purpose of reducing the size of the following expression, We shall 
denote 

E = e-'P% D= S = 4^^, S^’' = S-f^S-f’'S (56) 

— m p — m 

We have 


G^{x,a) 


(gauge covariant part, see eq. (21)) + 




(57) 


G^{x,^ 


= (gauge covariant part, see eq. (26)) + 

-i{A„W,A„(x)}^g^(£S) + 

1 f / d^E d^S 

^^d^duAp{x) yoppdp^dpp^ ^ ^dppdp^dpp 


+ 


+ ^{a,A.(a;),A,(x)}jr^(E^[H) + 

A / ^ A / M /■ /^ 1 ^ 1 d^E dS 

+ [5/iAy(x), Ap(a:)] / I -- - - S + -- - - 

^ Jp\2Adppdpi,dpp Sdppdp^dpp 


-Le. 


12 dppdp^dpr 


+ 


+ X2dpAi,{x)Ap{x) 
+ X3Ap{x)dpA^{x) 


[A 

Ip dpp 

f d 


Ip dpp 


(^ES’^p'^ 

ES^A 


+ 


+ 
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1 d^E 


+ K^{x)K{x)Kp{x) 


6 dppdp^dpp 


S- 


1 d‘^E dS 

8 dppdp^ dpp 


1 d^E dS 1 d^E dS 


4 dppdpp dpi 
1 d 


2dpp 


ES'^p 


8 dpudpp dpp 

1 d 


2dpp 


Est^^ 


^-E 

3 dppdpudpp 

+ ... 


( 58 ) 


While not all of these formulae look explicitly like boundary terms, perform¬ 
ing the integrations by parts shows that they are indeed. The boundary 
terms that contain higher number of fields A and their derivatives can be 
omitted in 4 dimensions, as they are either boundary terms of convergent 
integrals or vanish for a = 0. 
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